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Abstract 

We study the volatility of the MIB30-stock-index high-frequency data from Novem- 
ber 28, 1994 through September 15, 1995. Our aim is to empirically characterize 
the volatility random walk in the framework of continuous-time finance. To this 
end, we compute the index volatility by means of the log-return standard devia- 
tion. We choose an hourly time window in order to investigate intraday properties 
of volatility. A periodic component is found for the hourly time window, in agree- 
ment with previous observations. Fluctuations are studied by means of detrended 
fluctuation analysis, and we detect long-range correlations. Volatility values are 
log-stable distributed. We discuss the implications of these results for stochastic 
volatility modelling. 
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1 Introduction 



In this paper, we study the stochastic properties of MIB30-stock-index volatil- 
ity. The Black-Scholes/Merton (BS/M) framework for option pricing is based 
on the assumption of constant volatility [1]. In practice, volatility is time de- 
pendent and the characterization of its temporal evolution is one of the main 
task of the increasing number of physicists working in the financial field [2]. 
This problem has been thoroughly studied both by economists [3] and math- 
ematicians [4]. However, the emphasis has been mainly given to analytically 
tractable problems. 

Here, we use different approaches in the attempt to empirically characterize 
the MIB30-stock-index volatility. Among the various possible methods, we 
focus on the stochastic continuous-time volatility approach to option pricing. 

The paper is divided as follows: in section 2, we present an outline of the 
stochastic-volatility theory; in section 3, the main empirical results are dis- 
cussed; finally, conclusions are drawn in section 4. 



2 Theory 



In continuous-time finance, stochastic volatility, a t , can be modeled by means 
of two stochastic differential equations [4], a two factor model: 

dP(t) = p(P,t)P(t)dt + atPtydw^t) (1) 



da t = a{a t , t)dt + (3(a t , t)dw 2 (t) (2) 

where p(P, t), a(at, t) and /?(cr t , t) are deterministic functions of the spot price, 
P(t), or of the stochastic volatility, a t , and of time, t; the volatility, a t , is a 
stochastic variable, W\ and w 2 are standard one-dimensional Brownian mo- 
tions with correlation d{w 1 w 2 ) = pdt for some constant p. The two processes 
are independent if and only if p = [4]. The problem is then to find a unique 
solution {P,cr) for the system of stochastic differential equations (1) and (2). 
In the original BS/M model, a and (5 vanish and p and a are constant. 

If volatility is a stochastic process, continuous riskless hedging in the sense 
of BS/M (using an option and the underlying asset) is not possible [5]. This 
claim is based on a theorem concerning multi-factor stochastic models. In the 
particular case of the two-factor model of eqs. (1) and (2), in order to form 
a continuous riskless hedge, a financial instrument with price fully correlated 
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to volatility would be necessary [5]. A clear and exhaustive introduction to 
multi-factor models can be found in Marco Avellaneda's tutorials [6]. 



If the ansatz of eq. (2) is accurate, the empirical analysis of stochastic volatil- 
ity should lead to the determination of the coefficients a(a t , i) and j3(cr t , t) [7], 
thus completely specifying its stochastic dynamics. In practice, this task can 
be very difficult, due to data incompleteness and to possible intrinsic math- 
ematical difficulties. For instance, more than one set of the coefficients could 
well reproduce the known statistical properties of the volatility time series. 



3 Empirical Study 

We have analyzed MIB30 high frequency data from 28 November 1994 up to 
15 September 1995. MIB30 is an official index of the Italian Stock Exchange, 
it is composed by the 30 Italian shares with the highest capitalization and 
trading volumes, and is recorded every minute. The data set is composed by 
over 80,000 data: 420 data for every trading day. Considering the series of 
index values, Pj T , where r = 1 min and j = 0, . . . , 83579, we divide our data 
into non overlapping intervals or time windows of length T. We choose a time 
horizon or time scale, At, which is an integer multiple of r. We compute the 
logarithmic returns related to every interval T as follows: 



where N is such that T = NAt. 

Financial practitioners define historical volatility as the standard deviation 
of the logarithmic returns [8]. Following them, we estimate the volatility for 
every time window as follows: 



r nAt 



log 




, n = 0,... ,N-1; 



(3) 
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where r is the mean value given by: 




^ N-l 



(5) 



If At is measured as a fraction of year, we can define the annualized volatility: 



In order to investigate intraday properties of volatility, we chose a minutely 
time horizon and an hourly time window. The results of the 1393 annualized 
volatility estimates are plotted in Fig. 1. 
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Fig. 1. Annualized volatility estimate; the time window is 1 h; the time horizon is 
1 rain. 



In Fig. 2 we present the power-spectrum-density estimate computed by means 
of the correlogram method [9]. The peak at f / f s = 1/7 is due to a daily 
periodicity of the volatility values. In fact, in a day there are seven trading 
hours and with an hourly time window we get seven volatility estimates per 
day. Indeed, intra-day volatility is U-shaped: it is higher at the opening and at 
the closure of the market. This fact has already been observed by economists, 
and it is also known in the physics literature [10]; it probably reflects the 
lower trading activity around noon. As a further remark, we are not able to 
detect any low frequency seasonality or observe a clear flicker behaviour at 
low frequencies [11], as our time series is less than one year long. 



In Fig. 3, an estimate of the volatility probability density function is given. 
We compare the experimental histogram with a log-stable distribution whose 
parameters, 8 and 7 are obtained from empirical data. Stable or Pareto-Levy 
distributions [12] have been introduced in the sixties in finance and economics 
[13] and their scaling properties have been recently investigated in relation to 
the S&P500 stock index[14]. Zero mean stable distributions are described by 
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Fig. 2. Power spectrum density estimate. The sampling frequency, f s , is equal to 
\ I hour. 

the following equation: 



where 5 = 1 and 5 = 2 give, respectively, the well known Cauchy and Gauss 
stable distributions. A random variable is said to be log-stable distributed if 
its logarithm follows a stable distribution. 

In order to determine the experimental points, the volatility range has been 
divided into fifty equal intervals (bins). From Fig. 1, it can be seen that there 
are two outliers. If the outliers are taken into account, a log-Levy distribution 
with exponents 5 = 1.6 and 7 = 0.13 gives an acceptable fit of the experimental 
data (x 2 — 75 with 47 degrees of freedom), whereas the log-normal fit with 
mean value and standard deviation drawn from data does not agree in the 
tail region (x 2 ~ 10 4 ). Conversely, if the two outliers are rejected as bad data 
points, and a new histogram is computed using twenty bins, the log-normal 
fit is as good as the log-Levy fit (for the log-Levy fit, x 2 is 29.2, whereas for 
the log-normal, it is 28.7; this time there are 17 degrees of freedom). 

In Fig. 4, the results of detrended fluctuation analysis (DFA) are presented. 
DFA is used to investigate the presence of correlations in time series [10], [11], 
[15], [16]. A volatility walk is defined by means of the displacement y(t): 
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Fig. 3. Probability density function estimate; open circles: experimental histogram; 
solid line: log-Levy fit with 5 = 1.6 and 7 = 0.13; dash-dotted line: log-normal fit 
with mean value = 2.34 and standard deviation = 0.41. 

V(0 = £>(<)■ (8) 

i=i 

The mean square fluctuation, F(t), around the average displacement is given 
by: 

F(t) = J(Ay(t)*) - (Ay(0) 2 , (9) 

where Ay(t) = y(t + t) — y(t ), and (•) is the average over all the initial steps 
t . 

F(t) follows the scaling law: 

F{t) oc t z . (10) 
If long-range correlations are present [17], the scaling exponent attains values 

The points in Fig. 4, top, have been computed according to Eq. (9) with 
overlapping subseries. Forty Fit) estimates are plotted with their error bars. 
Errors are computed assuming a Gaussian distribution. The solid line is a 
least-square linear fit of the first six points, corresponding to one trading day; 
its slope gives z = 0.64 ± 0.23. The dash-dotted line is a least-square fit of 
the next thirty-four points and has a slope z = 0.76 ± 0.16. Therefore, from 
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Fig. 4. DFA; top: F(t) for the volatility walk; bottom: as top, but for de-U-shaped 
volatility walk, see text for explanation. 

our data we cannot safely conclude that there is an exponent cross-over at 
t — 1 trading day, as was found in ref. [11]. However, it is possible to argue 
that z > 1/2; indeed, a linear fit of the forty points gives z = 0.72 ± 0.10. 

In Fig. 4, bottom, we present F(t) computed by "de-U-shaping" volatility 
values according to the following recipe: 

a*(i)=a(i)/n(i) i = 1, . . . , 1393; (11) 

where the normalization coefficient, n(i), is given by: 

„ w = i sK^+^ i - (12) 
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where q = 1393/7 and n(i) has period 7: 

n(i + 7h) = n{i) h e Z. 



(13) 



In this way, the daily periodicity is removed and the intra-day estimate of z 
is 0.72 ±0.26, whereas the extra-day estimate is z = 0.77 ±0.16. The possible 
cross-over seems to be suppressed, and the full linear least-square fit gives 
z = 0.76 ± 0.10, that is z > 1/2. 



4 Summary and Conclusions 

In this paper, we have empirically studied the intra-day statistical proper- 
ties of stochastic volatility for the MIB30 index. Stochastic volatility poses 
serious problems in contingent claim analysis and risk management. In many 
cases, risk managers consider a daily time window for computing volatility and 
perform daily hedging. However, the preliminary analysis of a recent liquidity 
crisis of a hedge fund seems to suggest that frequent intra-daily hedging could 
be necessary [18]. 

Empirical analysis constrains the form of the stochastic differential equations 
describing the time evolution of volatility. 

In particular, we find in the MIB30 index that the volatility has a periodic be- 
haviour with a one trading day period. Due to the limited amount of available 
data, we are not able to detect any low-frequency seasonality. 

Probability-density-function estimates indicate that volatility is log-stable 
distributed; if outliers are taken into account a log-Levy distribution gives a 
better fit than the usually assumed log-normal distribution [7]. 

Finally, the DFA results are compatible with the presence of long-range volatil- 
ity correlations, but we cannot safely conclude that there is a crossover between 
intra-day and extra-day scaling exponents. 
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